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ABSTRACT 

We present a semi-analytical hydrodynamical model for the structure of reconfine- 
ment shocks formed in astrophysical relativistic jets interacting with external medium. 
We take into account exact conservation laws, both across the shock front and in the 
zone of the shocked matter, and exact angular relations. Our results confirm a good 
accuracy of the approximate formulae derived by Komissarov & Falle (1997). However, 
including the transverse pressure gradient in the shocked jet, we predict an absolute 
size of the shock to be about about twice larger. We calculate the efficiency of the 
kinetic energy dissipation in the shock and show a strong dependence on both the 
bulk Lorentz factor and opening angle of the jet. 

Key words: galaxies: jets - shock waves. 



1 INTRODUCTION 

Geometry (cross-sectional size, opening angle, and bend- 
ing) of supersonical, light jets is regulated, in general, by 
a complex system of oblique shocks. At certain circum- 
stanc es they take form of reconfinement shocks (|Sanderd 
1983). Such shocks have been considered to be responsi- 
ble for non-thermal activity in AGN radio cores (see, e.g., 
iDalv fc Ma rschcr 1988; Komi ssarov fc Fallelll997l - hereafter 



KF97; IStawarz et al\\200ai ) and , on much larger distances , 
in kiloparsec-scale radio knots (|Komissarov fc Falldfl998h . 
They are also predicted to operate in ma ssive X-ray binary 
systems llPerucho fc Bosch-Ramonl I2008T ) and in GRB col- 
lapsars (|Bromberg fc Levinsonll2007l V 

A direct way to verify whether a given source can be 
interpreted in terms of a reconfinement shock is to deter- 
mine whether location and extension of the source is con- 
sistent with a power of a jet and pressure/density of exter- 
nal medium. Under several simplifying assumptions a nalyt- 
ical fo rmulae relating these quantities were derived bv lFallei 
|l99ll ) and, for relativistic shocks, by KF97. We have devel- 
oped a semi-analytical model based on exact conservation 
laws and an exact dependence of a shock structure on an 
initial opening angle of a jet. Like in KF97, we adopt the 
cold jet approximation, i.e. we neglect the internal energy 
of the unshocked jet matter. The aim of this paper is to test 
the accuracy of the analytical formulae and to study the 
effects of including a transverse pressure gradients in the 
post-shock zone. 



Our models are described in §2. They are compared with 
analytical results of KF97 in §3. We study the efficiency of 
energy dissipation in the reconfinement shocks in §4 and 
discuss their possible 'astrophysical appearance' in §5. Our 
main results are summarized in §6. 



2 DESCRIPTION OF THE RECONFINEMENT 
MODELS 

The models we develop here are stationary, axisymmetric 
and purely hydrodynamical. We use a cylindrical coordi- 
nate system originating at the central source, with z-axis 
aligned with the jet symmetry axis and r denoting the cylin- 
drical radius. At every point, the flow is characterise d by the 
following parameters: bulk Lorentz factor F = 1/ yl — (P, 
rest-density p, pressure p and the angle between the velocity 
vector and the z-axis 8. We use the equation of state for the 
ideal gas 



(7-l)e, 



(1) 



where e is the internal energy density and 7 is the nu- 
merical coefficient, which for non-relativistic and for ultra- 
relativistic gases coincides with the adiabatic index with 
a value 5/3 and 4/3, respectively. For intermediate cases 
and/or mixtures of non-relativistic and relativistic gases 7 
takes intermediate values which might somewhat differ from 
the respective values of the adiabatic indices (see KF97) . 

In a stationary flow conservation of mass, energy and 
momentum is expressed by the following equations: 
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di{pu) = 0, 



(2) 
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Figure 1. Structure of the reconfincment shock for static external 
medium. 



diT'" = ; 
where 

1 — wu u + pg 

is the energy-momentum tensor, 



w — pc 



e+p 



(3) 



(4) 



(•>) 



is the enthalpy, u M is the 4-velocity (its spacial components 
are u 1 = F/3 1 ), and g^ v is the metric tensor of signature 

( (- ++). For the purpose of this study we assume a flat 

Minkowski space. 

The jet is modelled as a spherically symmetric adiabatic 
outflow from the central source into the cone of half-opening 
angle Oj. Equations ([2] — [3j) can be used to show that mass 
and energy fluxes of the upstream flow through any given 
solid angle must be conserved: 



d 



— [R 2 (r jW]Uj )] 



0; 



(6) 
(7) 



where R — \/r 2 + z 2 is the radial distance from the central 
source. The quantity FjWj/pj is invariant along R. 

In the cold jet approximation pj is negligible, so wj ~ 
PjC 2 . Then Fj is invariant along R, so from equation (0 we 



The total power of a jet is 
cos QAR 2 . 



(8) 



parameters for every point within the jet. 

In the interaction between the a jet and external matter 
a double oblique shock structure forms, but when external 
medium is static, it degenerates into a single shocked zone 
(see Fig. [TJ . In this scheme, the jet is bounded by the inner 
shock surface r s (z) and the shocked gas is bounded by the 
contact discontinuity r c (z). We denote the inclination angles 
of these surfaces by tana s ( c ) = dr s r c \/dz, respectively. The 
flow parameters (F, p, p and 9) are marked with the following 
subscripts: j - for the jet matter at r s , s - for the shocked 
matter at r s , c - for the shocked matter at r c and e - for 
external medium at r c . 

The conservation laws must be satisfied across the shock 
surface. Eq uations ([2] -EH) can be use d to obtain shock jump 
conditions (|Landau fc Lifshitd[l959h : 

M =0, (9) 

[pux] = , (10) 

[toul+p] = 0, (11) 

[Twu ± ] = ; (12) 

where mm and u± are the tangent and normal components 
(with respect to the shock surface) of the local velocity field, 



respectively. At the shock front r 3 , they give: 

P s cos(# 3 — a a ) = Pj cos(9j — ce s ) , (13) 

u 3 p s sin(# 3 — a s ) = UjPj sin(#j — a s ) , (14) 

u 2 s w a sin 2 (# s ~ a s ) +p a = u 2 Wj sin 2 (9j — a a ) + pj , (15) 

TsUsWs s'm(9 a — ce s ) = TjUjWj sin(#j — q s ) . (16) 

For the contact discontinuity there must be no flow through 

the surface, so the constraints derived from equations ((2]- 
[31) are much more simple: 



Pe ■ 



(17) 
(18) 



The purpose of our models is to find the geometrical struc- 
ture and the physical conditions of the shocked zone, given 
the conditions in the jet and in the external medium. 



2.1 Model 1 

In our first model we ado pt an assumption made by 
l|Bromberg fc Levinsonl |2007j) , that the shocked zone has 
no transverse structure, which means that for a given z: 
T c = T s , p c = p s , p c = p s , 9 C = 9 S . Knowing that p a = p e , 
we can solve equations 13 1 — 1 16 p for the 4 unknowns: F s , p s , 
9 3 , a s . This can be done explicitly using exact analytical for- 
mulae (see Appendix|X]|- We may then find r s by numerical 
integration over z. 

We considered also finding the contact discontinuity 
surface r c , by noting that a c = 9 S . But when we calcu- 
lated the mass flux across the shocked zone, we found that 
it is not conserved. Moreover, when neglecting the trans- 
verse pressure gradients, one cannot satisfy the transverse 
momentum balance needed to account for the curvature of 
the streamlines. 



2.2 Model 2 

A simple transverse structure of the shocked zone can be 
included in our model by assuming that parameters at op- 
posite boundaries are independent. For a given z we now 
have 4 more unknown parameters: F c , p c , a c and p s . There- 
fore we need 4 additional equations that would tie the flow 
parameters at the shock surface to the flow parameters at 
the contact discontinuity. 

We use conservation laws across the shocked zone in- 
troduced bv lBromberg fc Levinsonl (|2007l ): 



(/ 

dz 

d_ 

dz 



up cos 9r dr 



+ p a u s n a k- 



= 0. 



T^r dr 



(19) 



= 0; (20) 



where n s and n c are the vectors normal to the shock sur- 
face and contact discontinuity surface, respectively, oriented 
outwards the shocked zone. Equation (|19|l describes conser- 
vation of mass, while the equation set (|20[) includes con- 
servation of energy (p = 0) and of two momentum com- 
ponents (radial - p, = r, longitudinal - p — z). Note that 
Tl k n c k = Pcni 7^ 0. We can solve these equations by reduc- 
ing them to a system of linear ODEs (see Appendix [Bj. 
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3 GEOMETRICAL PROPERTIES OF 
RECONFINEMENT SHOCKS 

The crucial characteristic of the reconfinement shocks is 
their length scale, which may be estimated observationally. 
KF97 provided simple analytic formulae, in which they con- 
nect geometrical properties of the shock surface to physical 
parameters, such as the external pressure p e , the total power 
of a jet Lj, and its bulk Lorentz factor Fj. They assumed 
that: the pre-shock plasma is cold (pj <C p<?)\ pressure be- 
hind the shock is equal to the external pressure (p s = p e ); 
the half-opening angle Oj is small; and the pressure balance 
at the shock front, given by equation (|15p . can be approxi- 
mated by the formula: 



p s = fiUj pjC 2 sin 2 (8j 



(21) 



with n = 17/24. 

Below, we rewrite their results, using slightly different 
notation. We launch the jet from the distance z = Zo- Let the 
external pressure profile be p e (z) = po(z / zo)~ v (we expect 
T) > 0). The shock surface should satisfy a boundary con- 
dition r s (zo) = zotanGj, where Qj is the jet half-opening 
angle. Then the shock surface equation is: 



r s (z) = 



fo 

5A 



zo) 



where 5 = 1 — rj/2, and 
A = 



TTPqC 



(22) 



(23) 



is a characteristic length scal^B- The reconfinement is found 
for r) < 2(1 + «o/A) at 

z r = zo [1 + 5— J . (24) 
The maximum width of unshocked jet, 



4 



A \(l + 5)z 
is achieved at 
— Zr 

Zm ~ (l + S) 1 / 6 ' 

The aspect ratio of the jet is given by 



Zr (i+sy+ s J ' 

The shock surface inclination at zo is 



e = e 



(-?) 



(25) 



(26) 



(27) 



(28) 



The half-closing angle (equal to minus the shock inclination 

at z r ) is 



© r = e 



(29) 



1 Note that in the Eq. ((23J we have A <x L 1 ^ 2 ^ 2 , while mod- 

~ \/2 1/2 

els of nonrelativistic jets predict A a L/ (3j (see Eq. (1) in 
iKomissarovl l ll994l) and refs. therein). The reason for the differ- 
ence is that the term Lj includes the flux of the rest energy, Mj c 2 , 
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Figure 2. The ratio of reconfinement position z r to the char- 
acteristic length A as a function of Results for Model 1 
(dashed lines) and Model 2 (solid lines) are shown for different 
equations of state of the shocked matter: -y s = 4/3 (black lines) 
and 7 S =5/3 (grey lines). 
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Figure 3. The aspect ratio of unshocked jet r m /z r , divided by 
the value 0j/4 predicted by KF97, as a function of FjQj. The 
lincstyles are the same as in Fig. [2] 



For the case of uniform external pressure (77 = 0, 
8 — 1) and a jet originating close to the central source 
(zo A) we find the shock to be parabolic, with very 
simple characteristics: z r ~ A, r m /z r — 6j/4, z m — z r /2, 
0o — Or — Qj- We have tested these relations in our two 
models, the results are shown on Figs. [2] [5j as a function 
of half-opening angle Qj for a fixed Tj = 10. Other param- 

s _1 , po = 10~ 2 dyn and 



eters used were Lj = 10 46 erg 



zo = 10 cm. The characteristic length for these parame- 
ters is A = 2.74 • 10 18 cm = 0.89 pc. 

A very good agreement between the results of Model 
1 and the analytical formulae results from the same value 
of the pressure behind the shock (p s = p e ). Deviations for 
TjOj > 1 reflect the small angle approximation employed in 
analytical formulae. Small but systematic deviations of z r 
from A result from approximate pressure balance equation. 

In Model 2 the pressure behind the shock is system- 
atically lower than p e , but it cannot be fitted to a single 
power-law function of z. This results in longer reconfinement 
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Figure 4. The ratio of the maximum jet width position z m to the 
reconfinement position z r as a function of The linestyles 

are the same as in Fig. [2] 



Figure 6. Dissipation efficiency e^iss as a function of FjQj, cal- 
culated for Fj = 10. The linestyles are the same as in Fig. [2] 



1.4 
1.3 
1.2 
1.1 
1 

0.9 - 



Model 1 (y s =4/3) 
Model 1 (y s =5/3) 
Model 2 (y s =4/3) 
Model 2 (y s =5/3) 



100% 



10% 



1% 



0.1% 



0.01% 



0.1 



r j=5 
r ]= io 



=20 
=40 



10 



0.1 



Figure 5. The ratio of the half-closing 
opening angle 8j as a function of Tj&j. 
same as in Fig. [2] 



angle r to the half- 
The linestyles are the 



Figure 7. Dissipation efficiency e^iss as a function of FjQj, cal- 
culated for Model 2 with -y s =4/3. Line colour indicates the value 



structures (by a factor of about 2.2). We have found that for 
small and intermediate half-opening angles: r m /z r < 0j/4, 
z m < z r /2 and, accordingly, Gv < Oj. For large angles the 
deviations from analytical predictions are more pronounced. 
Nevertheless, the effects of independent values for the p s are 
not particularly strong. The analytical formulae are still very 
useful within the order of magnitude accuracy. 



4 ENERGY DISSIPATION 

The kinetic energy flux through the shock front is dissipated 
with efficiency 



(32) 



^kin(j) 

where 

F kin = pc 2 u ± (r - 1) 



(30) 



(31) 



and «i is the 4-velocity component normal to the shock 
front. Combining Eqs. (O, (J3TJ and {T4} gives 



As averaged over the entire shock front area, the efficiency 
of energy dissipation is found to strongly depend on the 



product Fj Qj . Results for both models with a fixed Tj 



10 



are shown in Fig. [6] We find that the averaged efficiency is 
very similar in both models and is insensitive to the value 
of 7 S . It approximately scales like ediss ~ 0.06(rjOj) 2 for 
FjOj < 1, but its increase slows down at TjOj > 1. 

In order to determine whether ediss is truly a function of 
VjQj, in Fig. [7]we present the results for Model 2 with 7 S = 
4/3 and different values of Tj. We find little discrepancy 
between the curves, which implies that it is a well defined 
dependence. 

We have investigated the ^-profiles of the dissipated en- 
ergy flux. In Fig. [8] we show the results for both models, 
with Tj — 10 and Qj = 5°. Although the reconfinement 
position z r is more than twice large in Model 2, as com- 
pared to Model 1, the total amount of dissipated energy is 
very similar. The dissipated energy profiles have a well de- 
fined maximum, which we denote as Zdiss,max- The ratio of 
Zdiss,max to z r is shown in Fig. [9] It is larger in Model 1, but 
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5 ASTROPHYSICAL APPEARANCE 
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Figure 8. Profiles of dissipated energy flux produced at the shock 
surface, calculated for Fj = 10 and @j = 5°. The linestyles are 
the same as in Fig. [5] 
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Figure 9. Ratio of the location of maximum of dissipated energy 
Zdiss max to the reconfinement position z r as a function of Tj@j. 
The linestyles are the same as in Fig. [2] 



larger than 1/2 in both models for FjOj < 1. It decreases 
strongly with increasing TjOj, for FjQj > 1. 

Noticing that the energy flux F w = wFu± = (pc 2 + 
'ye)Fu± is conserved across the shock front (see Eq. I16[) . 
one can find that, for jj ~ 7s, the efficiency of the internal 
energy production is 



F, 



e(s) 



F, 



Fh 



1 

7» 



(33) 

kin(j) 

A fraction of this energy is tapped by particles accelerated 
to relativistic energies and lost by nonthermal radiation. 
Such processes, if efficient, may significantly affect the shock 
structure. 

It should be noted that in the case of particle accel- 
eration with a broad energy distribution, most relativistic 
electrons may lose energy very efficiently even if the average 
energy dissipation efficiency is low. It means that, indepen- 
dently of the total energetics, the emissivity of such electrons 
will be maximized very close to the shock front and its spa- 
tial distribution will match the distribution of the energy 
dissipation. 



As theoretical analyses and numerical simulations demon- 
strate, formation of reconfineme nt shocks is ac companied by 
formation of reflection shocks l|Sanderslll983l ; KF97). Fur- 
thermore, depending on a distribution of pressure or den- 
sity of external medium, reconfinement and reflected shocks 
can form more or less abundant sequences of reconfinement 
shocks. Their radiative appearance is commonly modeled by 
assuming proportional ity of the emissivity to the gas pres- 
sure ( e. a. lG6meal2002j ; KF97). This leads to the predictions 
that most of the nonthermal radiation is produced around 
the reflection shocks. However, proportionality of the emis- 
sivity to the pressure is not what should be expected, if the 
efficiency of particle acceleration scales with the efficiency 
of energy dissipation. The latter is maximized at the shock 
fronts and, therefore, radiation emitted by most relativis- 
tic electrons will match geometrical structure of the shock 
fronts rather than the volume distribution of the pressure 
in the post shock flows. Of course, 'the shock front radia- 
tion' is likely to be accompanied by emission from the entire 
post-shock volume, by both slowly cooling lower energy elec- 
trons and by electrons accelerated in turbulent plasma in the 
2 nd order Fermi process. Specific geometrical and kinemat- 
ical structures of reconfinement shocks are expected to be 
reflected in polarization properties, provided that magnetic 
fields are dominated by the shock comp ressed random field 
(|Laing||l98dl : ICawthorne fc Cobblll990r ). This may explain 
perpendicular to the jet direction of pola rization (EVPA) o f 
radio knots in AGN kiloparsec scale jets (jBridle et qZ.lll994 V 



6 CONCLUSIONS 

- Semi-analytical models were developed to calculate the 
structure of reconfinement shocks based on exact conser- 
vation laws and exact angular relations. The approximate 
analytical formulae of KF97, that describe a shape of the re- 
confinement shock and its dependence on the power of a jet 
and the pressure of external medium, were confirmed with 
a very good accuracy, even for FjOj up to a few. However, 
the absolute size of the structure is found to be larger by 
a factor about two, when including the transverse pressure 
gradient in the post-shock flow. 

- The efficiency of energy dissipation in the relativistic 
reconfinement shocks scales approximately as (r^O^) 2 for 
TjOj < 1 and reaches about 6% at TjQj = 1. For both 
models, with or without the transversal pressure gradients, 
the efficiency is very similar and for a given value of TjOj 
practically does not depend on the bulk Lorentz factor. 
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APPENDIX A: SOLVING THE SHOCK JUMP 
EQUATIONS 

We show here a method by which the parameters of matter 
behind the shock may be determined in an exact analyti- 
cal manner from equations (|13l - [T6"Tl . First, we express the 
angular parameters with non-angular ones. From equation 
(|13[) we find an expression for a s , which is also a differential 
equation for the shock surface: 



dr s 



f3 s cos ( 



13 j cos t 



tana s = = - - (Al) 

dz p 3 sin 9 3 — f3j sin 9j 

From equations (|14p and (|A1|) . we find the deflection angle 
of velocity field: 

cos(6» s -9j)= . H , ' ' (A2) 



Now we find two more equations for two unknown parame- 
ters: r s , p s : 

w 3 _ Tjp s 



T 3 pj 

(rj - r^pswj = (pj 



(A3) 



(A4) 



)r,(r.p, + r i p i ) 

Equation (|A3|) is the result of dividing equation (|16|) by 
equation (fl4|) . Equation (|A4|| is derived from equation (fl~5)) 
by eliminating trigonometric functions using equations (j 131 
- I14[) and eliminating w 3 using equation (|A3|) . Using the 
equation of state for the shocked matter (and choosing the 
value of 7s), we finally find from equations (|A3[) and (|A4|I a 
quadratic equation for F s : 



hsPeiWj - Pj +Pe)] r 2 s + 

+ [(7.-i)(pj -p e )p^ 2 ] r J r s + 
+ {-w 3 (( ls ~i) Pj + Pe )]r 2 



(A5) 



0. 



Analyzing the parameters of this equation, we know that 
there is always only one positive solution. During our calcu- 
lations, we set an alert for unphysical F s < 1, but it never 
triggered. Finding T 3 , we calculate p 3 from equation ()A4|) 



and then we find the angular parameters: 9 S from equation 
(|A2|) and a 3 from equation (|A1|) . 



APPENDIX B: SOLVING THE EQUATIONS 
FOR CONSERVATION LAWS ACROSS THE 
SHOCKED ZONE 

The normal vectors ft 3 and n c are given explicitly by: 

n s — — cos a 3 e r + sin a s e z , 
n c = cos a c e r — sina c e z . 

Equations (|19I - |2TT)) may be expanded into: 



(Bl) 
(B2) 



d 

dz 

d_ 

dz 

d_ 

dz 



d 

dz 



up cos 9r dr 



Tuw cos 9r dr 



(u s p 3 sin 5 3 



— (r 3 u 3 w 3 smS 3 ) 



cos a 3 



u w sin 9 cos 9r dr 
(u 2 w 3 sin S 3 sin 9 3 + p 3 cos a 3 

[u z 'w cos 2 9 + p) r dr = 
(u 2 w s sin 8 3 cos 9 3 — p 3 sin ct 3 



(B3) 
(B4) 

(B5) 



PeTc 



(B6) 



+ p e r c tana c . 



To perform the integrals we have to describe the flow 
parameters between r 3 and r c as functions of r. We notice 
that the expressions to be integrated are linear functions of 
p and p (since the enthalpy w is also their linear function) 
and non-linear functions of T and 9. We decompose the in- 
tegrated functions into f(r) = g(T(r), 9(r)) ■ h(r) ■ r, where 
h(r) is one of p(r), p(r) or w(r). We assume that h(r) is 
linear: 



h(r) = h 3 + 



h c — h 3 



(r - r 3 ). 



The integrals are approximated with 
' fdr- 9(Fs,9 3 )+g(r c ,9 c ) 



h(r)r dr. 



(B7) 



(B8) 



Substituting these formulae to equations (|B3l - [B6)l we ob- 
tain a system of 4 differential equations for 8 variables: r s , 
p s , 9 3 , a 3 , p s , F c , p c and a c . The system is then closed by 
including differential forms of equations (|13l - fl6)l . 



